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PREFACE 

This  report  constitutes  the  final  report  of  Air  Force  contract 
F40600-82-C-0005,  Study  of  Asymptotic  Theory  of  Transonic  Wind  Tunnel 
Wall  Interference.  This  effort  was  conducted  under  the  sponsorship  of 
Arnold  Engineering  Development  Center,  ( AEDC ) ,  Air  Force  Systems 
Command  (AFSC),  Arnold  Air  Force  Station,  Tennessee  37389.  Dr.  Keith 
Kushman  was  the  AEDC  technical  representative  for  the  contract.  The 
manuscript  was  sumbitted  for  publication  on  May  8,  1984.  The 
reproducibles  used  in  the  reproduction  of  this  report  were  supplied  by 
the  authors. 
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1.0  INTRODUCTION 

The  problem  of  obtaining  free-field  transonic  characteristics  of  wind 

tunnel  models  will  continue  to  be  of  central  importance  to  aeronautical 

technology  for  the  indefinite  future.  Although  the  wind  tunnel  interference 

problem  has  received  considerable  attention  at  subsonic  speeds,  and  a  number 

of  classical  theories  have  been  developed,  e.g„,  Refs.  1-3,  there  Is  room 

for  progress  even  in  this  linear  regime,  as  exemplified  by  the  recent  efforts 
4 

of  Kraft  ,  who  devised  an  ingenious  procedure  based  on  Cauchy's  integral 
formula  for  two-dimensional  flow  over  confined  airfoils.  This  method  has 
been  extended  to  the  three-dimensional  flow  by  Sickles  and  Kraft  in  Ref.  5. 
Among  its  advantages,  In  contrast  to  the  older  concepts,  is  that  by  measurement 
of  two  flow  variables  it  eliminates  the  need  for  both  knowledge  of  the  wall 
characteristics  and  analytical  synthesis  of  the  model.  Presently,  effort  is 
underway  to  extend  this  procedure  to  nonlinear  transonic  flows. 

In  connection  with  the  transonic  regime,  other  procedures  have  been 
developed  which  possess  attractive  features.  One,  which  Is  of  great  interest, 
is  a  scheme  that  has  been  developed  by  Murman.  This  "post  test  assessment 
method",  which  has  been  implemented  in  two  dimensions  in  Ref.  6,  is  similar 
in  some  respects  to  an  approach  developed  earlier  by  Kemp7,  and  assumes  a 
knowledge  of  the  experimental  pressures  on  the  model  and  walls.  It  uses 
modern  computational  and  optimization  procedures  to  determine  whether  the 
tunnel  Mach  number  and  model  angle  of  attack  are  correctable  in  the  sense 
that  almost  free-field  model  pressures  can  be  obtained  at  practically  altered 
values  of  these  parameters.  For  this  purpose,  an  inverse  problem  needs  to 
be  solved.  In  three  dimensions,  surface  pressures  over  the  model  are  generally 
not  available,  and  efforts  are  currently  underway  to  modify  the  method 
toward  the  use  of  less  information  regarding  the  model  near  field.  Examples 
of  schemes  of  this  kind  are  given  in  Refs.  5,  8,  and  9.  If  this  goal  can 
be  achieved,  this  process  will  be  of  value,  including  treatment  of  cases  in 
which  the  walls  are  relatively  close  to  the  model.  In  spite  of  the  potential 
utility  of  this  method,  there  is  a  need  for  approaches  that  can  reduce  the 
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number  of  Input  parameters  necessary  to  compute  the  correction,  shed  light  on 
the  physics  of  the  wall  interference  phenomena,  simplify  the  necessary  compu¬ 
tations,  and  be  generalized  to  three  dimensions,  as  well  as  unsteady  flows. 
Asymptotic  procedures  provide  such  advantages.  Furthermore,  they  can  provide 
valuable  interactions  with  the  other  methods  previously  mentioned  to  suggest 
possible  improvements  as  well  as  deriving  beneficial  features  from  than  as 
well.  For  example,  one  concept  presently  being  considered  in  the  simplifi¬ 
cation  of  the  Murman  method  of  Ref.  6  for  three-dimensional  applications  Is 
unfolding  a  singular  character  of  the  near  field  rather  than  obtaining  the 
model's  shape  from  an  inverse  solution  of  a  problem  involving  specified 
pressures  based  on  measured  values.  In  Ref.  10,  an  asymptotic  procedure  for 
two-dimensional  transonic  flow  was  developed.  From  this  analysis,  the  singular 
character  can  be  obtained  from  certain  limit  processes.  Moreover,  nonlinear 
integral  theorems  as  well  as  the  asymptotic  structure  of  nonlinear  integral 
equations  arising  in  the  matching  scheme  occurring  in  the  asymptotic  analysis 
could  be  of  use  in  the  procedure  of  Refs.  4  and  5. 

Other  potential  applications  of  the  asymptotic  theory  involve  adaptive 
wind  tunnels.  An  example  of  one  configuration  is  described  in  Ref.  11. 

The  analysis  given  herein  will  be  oriented  toward  the  three-dimensional 
generalization  of  asymptotic  solutions  developed  for  the  two-dimensional  case 
in  Ref.  10.  The  latter  represent  the  application  of  perturbation  theory  to 
the  solid  wall  Interference  problem  at  transonic  speeds.  Whereas  many  of  the 
previous  methods  can  handle  arbitrary  wall  to  model  dimension  ratios,  h,  the 
perturbation  procedure  assumes  h  to  be  large.  This  approximation  is  useful 
for  many  practical  cases  in  which  it  Is  desired  to  minimize  the  wall  inter¬ 
ference.  Furthermore,  even  for  situations  where  h  is  not  so  large,  the 
expansions  appear  to  have  extended  validity. 

A  previous  analysis  along  these  lines  was  conducted  by  Chan  ,  who 
treated  the  two-dimensional  porous  wall  case  based  on  asymptotic  developments 
similar  to  those  given  for  transonic  lifting  line  theory  in  Ref.  13.  Because 
of  Its  interest  in  .connection  with  compliant  boundary  applications  and  the 
fact  that  the  Chan  solutions  apparently  do  not  subsume  it,  we  have  analyzed 
the  solid  wall  case.  By  contrast  to  the  method  of  Ref.  12,  we  employ 
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"intermediate  variables"  to  match  the  model  near  and  far  fields.  Although 

14 

slightly  more  cumbersome  than  the  approach  utilized  by  Van  Dyke  ,  it  pro¬ 
vides  a  reliable  means  of  ensuring  that  all  the  proper  terms  are  included 
in  both  representations. 

In  addition,  numerical  solutions  have  been  obtained  for  this  theory, 
as  typified  in  Ref.  10,  to  give  some  insight  into  the  nature  and  magnitude  of 
the  interference  effects.  The  work  of  Ref.  12  in  this  sense  was  strictly 
formulational,  with  no  computational  application  or  quantification  of  the 
interference  given. 

This  report  provides  three-dimensional  generalizations  of  the  two- 
dimensional  theory  of  Ref.  10.  The  development  of  these  models  was  conducted 
under  Air  Force  Contract  F40600-82-C0005  and  monitored  at  Arnold  Engineering 
Development  Center.  The  program  consisted  of  the  following  two  tasks: 

Task  1 .  Asymptotic  procedures  will  be  considered  for  two  limiting 
cases:  slender  Configurations  representative  of  fighter  aircraft, 
and  high  aspect  ratio  configurations  representative  of  transports. 

The  feasibility  and  general  approach  of  applying  asymptotic  theory 
to  these  cases  will  be  determined. 

Task  2.  Based  on  the  results  of  Task  1  and  following  consultation 
with  the  sponsor  of  the  work,  either  the  slender  or  high  aspect  ratio 
case  will  be  selected  for  development  of  the  theory  for  assessment/ 
correction  of  wall  interference  for  three-dimensional  transonic  flow. 

On  the  basis  of  Task  1,  the  low  aspect  ratio  received  emphasis  in  this  progra'm. 
Within  this  framework,  the  case  of  a  slender  fighter  vehicle  confined  within 
cylindrical  wind  tunnel  walls  was  treated.  The  corresponding  theory  described 
herein  provides  the  formulation  of  a  numerical  problem  whose  solution*  gives 
the  wall-induced  interference  correction.  Although  the  results  apply  to 
sting-mounted  models  confined  by  cylindrical  solid  walls,  they  can  be  readily 
extended  to  other  support  arrangements  and  wall  cross  sections.  These 
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generalizations  simplify  the  numerical  work  required  for  the  assessment/correction 
as  compared  to  purely  computational  schemes.  Porous  and  slotted,  as  well  as 
other  boundary  conditions  such  as  a  specified  pressure  distribution  on  a  control 
surface  for  interference  assessment,  and  adaptive  applications  can  also  be 
handled.  The  theory  applies  to  the  situation  when  the  characteristic  lateral 
dimension  of  the  model  is  small  compared  to  its  length,  and  the  tunnel  height 
is  inversely  proportional  to  this  lateral  dimension.  To  be  provided  in  this 
report  is  information  from  the  theory  on  the  forces  and  pressures  associated 
with  the  interference. 

In  another  phase  of  the  effort,  progress  has  been  made  in  the  development 
of  a  comparable  model  for  confined  high  aspect  ratio  shapes.  Basic  ideas  have 
been  worked  out  which  represent  a  generalization  of  transonic  lifting  line 
theory.  Means  have  been  Identified  which  will  be  helpful  in  matching  the 
vortex  sheet  far  field  representation  consisting  of  a  lifting  line  reflected 
in  the  walls  with  the  nonlinear  almost  two-dimensional  near  field.  These 
developments  will  be  summarized  in  what  follows.  On  the  basis  of  this  effort, 
recommendations  for  future  study  will  be  made. 
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2.0  ANALYSES 

2.1  CONFINED  SLENDER  CONFIGURATIONS 


2.1.1  Overview  and  Roadmap 


In  what  follows*  a  slender  confined  airplane  model  shown  schematically 
in  Fig.  1  will  be  considered.  The  flow  structure  consists  of  three  zones 
in  which  different  approximations  for  the  perturbation  potential  apply. 

These  regions  are  indicated  in  Fig.  2.  Near  the  axis  of  symmetry  of  an 
equivalent  body  of  revolution  having  the  same  streamwlse  distribution  of 
cross-sectional  area  as  the  complete  airplane  {axis  layer),  lateral  gradients 
dominate.  In  a  "central  layer",  if  a,  the  angle  of  attack,  and  the  charac¬ 
teristic  thickness,  6,  are  such  that  a/6  =  0(1),  as  6  *►  0,  the  flow  is 
nearly  axi symmetric  and  can  be  characterized  as  a  nonlinear  line  source. 
Asymptotic  representations  for  the  central  and  axis  layers  exist  in  which 
the  first  order  terms  are  those  associated  with  the  unconfined  flow.  The 


Fig.  1  Slender  vehicle  confined  inside  cylindrical  wind 

tunnel  walls  indicating  Cartesian,  cylindrical,  and 
spherical  polar  coordinates  used  in  analysis. 
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Fig.  2  Front  view  of  wind  tunnel  model  confined  by  cylindrical 
walls,  showing  regions  of  applicability  of  asymptotic 
expansions. 

second  order  corrections  of  these  regions  are  due  to  the  wall  effects.  A  third 
region  denoted  as  the  wall  layer  can  be  identified,  where  the  assumption  of 
small  wall  perturbations  is  invalid.  Here,  other  simplifications  apply  which 
represent  the  slender  airplane  as  a  multipole  reflected  In  the  walls. 

The  effect  of  the  walls  on  the  flow  field  is  deduced  by  solving  the 
second  order  problem  for  the  central  layer.  This  consists  of  the  equation 
of  motion ,  hereinafter  referred  to  as  the  "variational  equation ",  subject 
to  boundary  conditions  devised  from  matching  the  wall  and  axis  layers. 
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Besides  furnishing  boundary  conditions,  the  matching  process  is  used  to 
determine  unknown  elements  appearing  In  each  of  the  representations  for  the 
various  layers.  It  consists  of  comparing  these  representations  in  a  common 
region  to  two  of  them. 

In  what  follows,  details  of  the  matching  process  will  be  provided. 

Once  the  asymptotic  form  of  the  solution  is  determined,  it  will  be  used  to 
evaluate  the  influence  of  the  walls  on  the  pressures  and  forces  on  the 
confined  body. 

2.1.2  Mathematical  Details 


Karman-Guderley  Approximation 

The  flow  over  a  slender  airplane  wind  tunnel  model  shown  in  Fig.  1  is 
considered  in  which  the  surface  of  the  test  article  is  given  by 

r  =  <5F(x,9)  ,  (1) 

in  cylindrical  coordinates.  In  the  notation  of  Fig.  1,  normalized  coordinates. 


(in  which  c  Is  a  characteristic  body  length)  are  introduced,  and  the  bars 
signify  dimensional  quantities. 

The  Karman-Guderley  (KG)  representation  of  the  velocity  potential  in 
small  disturbance  theory  is  given  by  the  following  asymptotic  expansion: 


=  x  +  <52<j>(x,r,0;K,H,A)  +  ,  (3) 

which  is  an  approximate  representation  of  *  in  the  KG  limit, 

1 

x,r  =  6r,6  ,  K  =  — —  ,  H  =  ,  A  =  j  fixed  ,  as  5  ♦  0  .  (4) 

61  c  0 
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Central  Layer  Expansion 

Within  the  KG  limit,  secondary  limits  involving  H  +  »  can  be  considered*. 
For  slender  airplane  models  confined  by  wind  tunnel  walls,  various  regions 
can  be  Identified  in  which  approximations  for  the  perturbation  potential  £  hold. 
These  approximations  are  expressed  as  asymptotic  expansions.  In  the  central 
region  depicted  in  Fig.  2  (away  from  the  x  axis  and  the  walls),  the  perturba¬ 
tion  potential  <J>  is  almost  axisymmetric  and  has  the  representation 

♦  =  *0(x,r)  +  (5a) 

which  holds  in  the  central  region  limit 

x,r  fixed  as  H  *  00  ,  (5b) 

where  the  quantity  y.  (H)^  is  a  "switchback"  term  inserted  for  matching. 

For  purposes  of  studying  the  possibility  of  obtaining  interference-free 
conditions,  the  KG  similarity  parameters  K  and  A  can  also  be  expressed  in 
the  perturbation  forms 

K  =  Kq  +  v-,  ( H)K]  +  ...  (5c) 

A  =  Aq  +  tc1(H)A1  +  ...  .  (5d) 

In  (5a),  (5c),  and  (5d),  the  flow  quantities  are  small  perturbations  about 
their  free  field  (subscript  0)  values  consistent  with  H  -►  «».  In  what  follows, 
only  (5c)  will  be  considered  to  obtain  values  of  K-|  for  interference-free 
loading  on  the  model . 


*The  approximation  scheme  contrasts  with  another  one  which  is  embedded  in  a 
full  potential  framework  under  development  by  the  authors. 
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On  substitution  into  the  KG  equations  and  retaining  like  order  terms, 

(5a)  and  (5c)  lead  to  the  following  equations  for  the  secondary  approximations 

* 

4>0>  and  t>-j : 

k-(rH)*0 1*0  +r(f*0.)  -0  <6a> 

L  xJ  xx  r  '  r/f; 

"  <v+i\x\  +  r  =  0  (6b) 

fKn-  (y+T  )+o  1*1  *  (y+1  )d>-|  *0  +r(^u)  ■  "Mo  ■  (6c) 

L  u  xJ  xx  x  xx  r  \  r/p  xx 

The  forcing  term  in  (6c)  is  necessary  to  achieve  interference-free  flow 
conditions.  Its  retention  requires  that 

^(Hj-v^H).  (7) 

For  matching,  a  procedure  described  in  Appendix  A  based  on  Green's  theorem 
gives  a  far  field  representation  for  4>q  for  a  "ray"  limit  in  which 

R  =  \/x2  +  Kq?z  -*■  »  ,  u  *  cos"1  ^  fixed  ,  (8a) 

where 

R  =  R//K^  =  \/  Xa  +  ra  ,  X  =  x//k^  .  (8b) 

From  this  process,  the  far  field  representation  is 
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A„  B, 


_9  _9  !2  /  2  \  A  /An  R\ 

4>n  ■  tj-  +  —  costa  +  —  cos3oi  +  |-  cosoj  +  —  P0(costa)  +  0[ — -)  , 
0  R  R2  r3  V  5  I  Rj  2  \  R**  / 


(9a) 


where  Pn(cosw)  is  a  Legendre  polynomial  and 


a  =  _  S(l) 
^  4tt 


,  4trBn  =  -S ( 1 )  +  f  S(x)dx  +  iLrtll  f  dxf  r<Ji£  (x,r)dr  , 
-  JQ  lX)  •'-»  JQ  ux 


c  sfnuxtn 

-2  IB^Kq 


(9b) 


*27 r 


S{x)  «  1-  P  F2de  . 
* 


(9c) 


The  constant  A  has  not  been  determined  and  will  not  require  an  evaluation  in 
what  follows. 

Wall  Layer  Expansion 

The  central  region  expansions  (5a),  (5c),  and  (5d)  are  nonuniformly  valid 
at  the  walls.  This  is  associated  with  the  idea  that  the  wall  Induced 
perturbations  are  no  longer  small  In  that  domain.  Instead,  the  appropriate 
representation  is: 

<j>  *  EQ(H)v>Q(x+,r+)  +  e-j (H)i^-j (x^,r"^)  +  •••  (10) 

which  is  valid  in  the  wall  layer  limit 

x+  =  ,  r+  =  ,  fixed  as  H  ®  . 
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By  the  substitution  procedure  previously  described,  the  wall  layer  approximation 
can  be  shown  to  satisfy  the  following  hierarchy; 


+ 7  (rtv)rt  ■  ° 

WxV  +  +  Wy.V 


(11a) 


(lib) 


and 


1=1,2,— 


(lie) 


Integral  Representation  for 

If  S C 1 )  ^0,  as  for  the  case  of  a  sting  support  rather  than  a  magnetic 

4. 

suspension*,  Vq  must  behave  as  a  reflected  source  near  the  origin  R  =0  in 
order  to  match  with  the  dominant  term  of  d>Q  given  in  (9a).  From  Appendix  B, 
the  appropriate  solution  with  this  property  can  be  derived  to  be 

^o'7=j^S9nX+tM|  •  (,2a) 

where 

Xf  =  x+//K^  02b) 

M=M0  +  M1  •  (12c) 


♦This  analysis  can  be  extended  to  treat  the  magnetic  suspension  case. 
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M 


0 


cos  kX+Kg(kr+)dk 


( 1 2d.) 


M 


K1(k)I0(kr+)| 

-  HrjikT-  j 


dk  , 


(12e) 


and  I  and  K  are  modified  Bessel  functions  of  order  n. 
n  n 

f 

Asymptotics  for  Small  R 

For  matching,  the  behavior  of  for  R^  =  R/H  +  0,  w  fixed,  is  needed. 

In  this  connection,  the  representations  ( 1 2d)  and  (12e)  can  be  expanded  near 
the  origin  by  Taylor's  theorem  or  another  method  involving  differentiation 
with  respect  to  r^  and  described  in  Appendix  C.  This  process  gives 


_  S(li 


— qr  +  an  +  bnR^  P, (cosco)  +  ••• 
41TR1*  U  U  L 


where* 


- — r 

TT2  J0 


k2 


Kj  (k) 


dk 


1  /•<»  k2K^  (k) 

4uz  Jq  ~W 


dk  . 


(13a) 


(13b) 


(13c) 


The  integrals  in  (13b)  and  (13c)  are  convergent,  can  be  evaluated  once  and 
for  all  numerically,  and  are  independent  of  XT  and  r  . 


*S1milar  quadratic  growth  to  the  R 


+2 


term  occurs  in  the  two-dimensional  case. 
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Hatching  of  Central  and  Mall  Layers 

An  intermediate  variable  and  limit  are  introduced  such  that 


Rn  =  nW  *  “  fixed  as  H  *  08  •  (14) 

Here,  n(H)  Is  of  an  order  between  unity  and  H.  This  property  is  expressed 
symbolically  as 

1  «  n{H}  «  H  .  (15) 

The  central  region  representation  for  $  in  the  intermediate  limit  (14)  is 


A0  Bq  Cq  1  2  )  A 

<j)  =  — * —  +  ■  ■  COSO)  +  •  *  {C0S3w  +  f  COSO)  /  + 

_ 2 n 2  -  3n  3  r  5  l 


nhn  n*R*  n3R •  ( 


_ _  •  p?(cosot) 

/ityi'R*  2 


©  © 


© 


(tnnRn , 

^rl  + VHlS  +  u1tHHl(Rn,“) 

©  ©  © 


(16) 


where 


A0  =  »  Bo  “  Vv%  *  C0  "  * 


The  central  and  wall  layer  representations  for  <j>  have  an  "overlap"  domain  of 
common  validity  as  can  be  demonstrated  by  writing  the  wall  layer  representation 
in  the  intermediate  limit  (14).  This  gives 
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* =  e°(H)  r  dr + a° + b°  £  W“->  +  -  ♦  ^  (& v) 


IT  n 


0  ©  © 


+  el(H>*’i(S'l,n-“) 


(17) 


where  an  additional  switchback  term  has  been  inserted  in  (10)  to 
match  .  Upon  identifying  common  terms,  denoted  below  by  «*,  matching 
demands  that  the  following  identities  hold: 


en  =  1/H 


A„  = 


1  sm 


©  *© 


e-|  (H)  =  1/H3 


*,  =  R+‘3 


I,  ?  A  4iTbn 

Cq(COs3(jJ  +  ^  cosw)  +  -j—  P2(costu)  >  -  a 

v^fT  I  vKa  “ 


as  R 


(18a) 

(18b) 


(19a) 


(19b) 


Equation  (19b)  is  consistent  with  the  fact  that  the  terms  on  the  right  hand 
side  are  solutions  of  (lib)  which  is  to  dominant  order  identical  in  form 
to  (6c).  These  terms  are  particular  and  homogeneous  solutions  identified 
in  (9a).  For  <0, ,  the  homogeneous  solution  can  be  obtained  from  differentiation 
with  respect  to  X  of  the  switchback  term<^  and  adjusting  the  quadrupole 
strength  for  matching.  The  quantity^  itself  is  obtained  from  differentiation 
of  vQ  with  respect  to  X.  It  satisfies  an  equation  identical  to  (11a)  and 
can  be  shown  to  be  -  BQ  costa  |r+"2  +  8tTbQRtj  as  R+  -*■  0. 
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y-,  (H)  =  1  /H 3  (20a) 

<J>t  =  bn  [MIL  R2P9(cosw)  +  STrBnR  cosw  -  — r— 1  as  R  -*■  00  .  (20b) 

1  2  u  /i^J 

In  the  formulation  of  the  numerical  problem  for  <j>.| ,  the  last  term  in  (20b) 
will  be  neglected  in  the  far  field  specification. 


(21a) 

(21b) 


noting  that  (6b)  is  identically  satisfied  by  the  constant  aQ. 


Axis  Layer  Expansion 

The  third  region  of  importance  is  a  zone  near  the  x  axis  shown  in 
Fig.  2  and  denoted  as  the  axis  layer.  In  this  region,  as  in  the  unconfined 
case,  cross  flow  gradients  dominate.  This  is  related  to  the  fact  that  the 
scale  of  the  r  gradients  in  this  domain  is  5.  The  asymptotic  representation 
for  the  velocity  potential  is 


jj  =  x  +  c[2(62*n«)^*(x)  +  S2<J>*(x,r*  ,9,K,H,A)  +  •■•]  (22a) 

which  is  valid  In  the  axis  limit 

r*  =  r/6 ,K,H,A  fixed  as  6  ^  0  .  (22b) 

Here,  r  and  9  are  cylindrical  coordinates  shown  in  Fig.  2. 
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As  a  representation  of  the  wall  effects,  the  perturbation  potential,  $*,  has 
the  development 

<J>*  =  +  ^(H)**  +  x1  (H)*^  +  —  (23a) 

which  holds  in  the  limit 

r*,K,A  fixed  as  S  ♦  0,  H  ®  independently  .  (23b) 

In  (23a),  the  quantity  Involving  is  a  switchback  term  inserted  to  match 
with  Its  counterpart  involving  in  (5a). 

Problems  for  <J>g,  and  <j>^ 

On  substitution  of  (23a)  into  exact  equations  and  boundary  conditions, 
the  boundary  value  problems  satisfied  by  $*,  and  <|>*  in  (23a)  can  be 
determined. 

For  **, 


i  a  ra» o\ .  i  a*»s  _ . 

r*  3r*  \  3r*  J  r*2  30a 


(24a) 


holds  in  the  region  external  to  the  cross  flow  boundary  B  =  0  shown  in  Fig.  3. 
This  boundary  is  the  projection  of  the  model  shown  in  Fig.  1  in  a  plane 
perpendicular  to  the  free  stream.  Equation  (24)  is  solved  subject  to  the 
flow  tangency  condition  on  the  boundary - 


an 


B 


FF. 


F2  +  Fg 


(24b) 


where  n  signifies  differentiation  in  the  normal  direction  to  the  surface 
B =  r*  -  F(x,9)  = 0,  and  B  signifies  evaluation  of  this  derivative  on  the  surface. 
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y* 


Fig.  3  Cross  flow  plane  geometry. 


On  the  basis  of  Gauss'  divergence  theorem  and  (24b), 


<j>*  -  jin  r*  as  r*  +  »  . 

(24) 

The  problems  for  the  other  quantities  are 

© 

ii 

tr 

<j 

ii 

$ 

(25a) 

%  _s*i  . 0 

3n  B  3n  B 

(25b) 

By  matching  considerations,  it  is  anticipated  that 

<  “  as  r*  -*•  . 

(25c) 
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Solutions  for  Axis  Layer  Approximations 

On  the  basis  of  Eqs.  (24)  and  (25),  the  following  multi pole  representations 

s*(x)  (01/  x  ^  g*{0)(x)cos  ne 

♦s-  a, )ll"'*+sr  ’(x)*!' 

(26a) 

tj>^  =  constant  = 

(26b) 

♦t  -  95°  >(*)  . 

(26c) 

where  without  loss  of  generalization,  lateral  symmetry  has  been  assumed  in 
the  sense  that  F{x,6)  =  F(x,-0). 

Matching  of  Axis  and  Central  Layers 

Following  a  procedure  similar  to  that  employed  in  the  subsection  relevant 
to  Central -Wall  Layer  matching,  a  suitable  Intermediate  variable  is  introduced. 
For  the  overlap  domain  of  common  validity  of  the  axis  and  central  layers,  this 
variable  is  assumed  to  be 


r  -  r 

rc  "  clfiT 


(27) 


where 


<5  «  c(«)  «  S'1  •  (28) 

The  axis  layer  representation  for  $  written  in  Intermediate  variables  is  then 
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J  *  x  +  2(62  Jln5l^*(x)  +  6! | [ln;r? - «]  £”*(x)  +  gg*0'(x) 

©  ©  ®  © 

+  Vi +  Tik(,)<x)  +  — 1}  • 

©  ©  ’ 


g>*  =  S'  (x)/2tt  . 


For  the  central  layer,  the  corresponding  representation  is 


jj  *  x  +  fia  J^x)j*n6  +  £ncrJ  +  gQ(x)  +  ^  g^x)  +  ••• 

©  ©  ©  ©  © 

On  the  basis  of  (29)  and  (30),  the  following  matchings  apply: 

©  +  ©»© 

&*m  -  •s'M 

©~© 

Already  matched. 


(29) 


(30) 


(31) 


Q"© 

gj$(0*(x)  =  gQ(x) 


(32) 
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CD"© 

T,.  -  1/H  •  (33a) 


©~© 
T,  =  1/H3 


(33b) 


(34) 


g 


* 

o 


(35) 


This  completes  the  matching  process. 

Suninary  of  Representations  for  Potentials  in  Various  Regions 

Collecting  the  previous  results,  the  asymptotic  expansions  for  the 
velocity  potential  In  the  Axis,  Central,  and  Wall  layers  are 

Axis  Layer 


4 

U 


x  +  l  {62Anfi)S'(x)  +  fi2 


*j(x,r*) 


0 

i r 


g^x) 

“h3- 


(36) 


r*  =  r/fi 


Central  Layer 


i 

u 


x  +  Sz|$0(x,r)  +  -ft-  +  -Jy  4>i  (x,r)  + 


r  =  fir 


(37) 
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Mall  Layer 


(38) 


/  =  *  P+  .  L 
H  ’  H  ’ 


...  s(i )  r(i  Ki(k) 
0  =  Jo  (k2  '  ^ 


dk 


Associated  with  obtaining  interference-free  loadings,  the  KG  similarity 
parameter  has  the  expansion 


K=K0+^K1  +  -  • 


(39) 


Interference  Pressure  Distribution  on  the  Model 

On  the  basis  of  Ref.  15,  the  pressure  distribution  Ca  on  the  model  is 


(S2 


£n<5)S"(x)  -  2<frJ 


+  +  < 


(40) 


where 


v 


* 

B 


w*  B 

WB  r^  30 


9 

r*=F 

B  " 

r*=F 


Substituting  (36)  in  (40),  this  pressure  can  be  written  as 
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Hence,  the  wall  interference  effect  on  the  surface  pressures  denoted  as 
ACp  is  independent  of  9  and  Is  given  by 

ACp  = "  *  (42} 

From  (41)  and  (42),  the  loading  on  the  model  is  constant  In  each  streamwise 
plane  x  =  constant.  Accordingly,  the  walls  induce  no  lift  perturbation  on 
the  model.  Only  a  drag  perturbation  occurs.  This  is  considered  in  the 
next  section. 

Tunnel -Wall  Interference  Drag 

The  drag  D  of  the  wind-tunnel  model  can  be  computed  by  applying  the 
analysis  of  Ref.  15,  which  was  performed  for  unconfined  flows.  Therein, 

D  is  computed  from  a  momentum  flux  in  the  x  direction  across  the  boundaries 
of  a  cylindrical  control  surface.  The  curved  boundaries  of  this  surface 
coincide  with  those  of  the  Axis  Layer.  From  this  procedure,  D  in  an 
unconfined  flow  is  given  as 

^■+  sHnsfWu)]  =  fij^Dgd)  -  \  jj^Vo »F,9)FFx(l ,e)de 


In  (43), 

g(x)  =  lim]<j>  -55'(x)4nr{  (44) 

r+O  (  ) 

q  =  dynamic  pressure  . 
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The  wall  Induced  perturbation  to  the  drag  AD  is  computed  from  the  change  in 
4>*  due  to  the  second  and  third  order  terms,  using  (43)  whose  structure  is 
otherwise  unaffected  by  the  walls.  The  result  is 

Sp-  f’s'UlgJtxJdx  .  (45) 

H  Hd  -'o  1 

Since  S  is  proportional  to  6  multiplied  by  the  unsealed  cross  sectional 
area,  the  change  in  drag  coefficient  normalized  to  the  product  of  the  model 
frontal  area  and  pp  ,  HACg  is  proportional  to  S2/H2.  This  quantity  is  in  turn 
proportional  to  the  blockage  ratio,  ar,  i.e.,  S2/h2 . 

Determination  of  Interference-Free  Conditions  for  Drag 

For  a  given  model  geometry  or  an  affinity  of  models,  Eq.  (45)  can  be 
written  in  the  form 

HACn 

f(KQ,K1,A,Y)  .  (46) 

“r 

This  relation  Is  a  similarity  law  In  which  H,ar  are  separated  out  of  the 
universal  variation  given  by  f.  The  latter  is  determined  solely  from  (45) 
through  g-j .  To  determine  this  quantity,  the  following  boundary  value 
problem  (PI)  needs  to  be  solved  for  the  previously  given  variational 
equation  for  4><| . 
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Pi: 

[k0  -  (y+1 ) <J>q  L  -  (y+IHq  + 

L  u  wx-l  XX  XX  X 

f  K)-  ‘  '  K’*°xx 

(6c) 

11m  =  0 

r+0  ? 

(47) 

$1  =  bgR^Uosod)  +  STrbgBgR  cosu  , 

R  ■+■  « 

(20b) 

b.  _  s(D  p  kZ|<i(k) 

(13c1) 

4irBn  =  -  S(l)  +  S(x)dx  + 

0  J0 

i*CO  /•  00 

I  dx  1  ?*£  (x,r)dr 

/_«  J  Q  ux 

(9b) 

where  c(>q  is  the  solution  of  (6a)  subject  to  (9a)  and 

11m  r<jjQ  =  S'(x)  for  0  <  x  <  1  and  <f>Q_  =  0  for  x  >  1  . 
r-0  ?  ? 

Once  the  problem  Pi  is  solved,  using  computational  methods  such  as  successive 
line  overrelaxation,  the  results  can  be  applied  to  extrapolate  to  zero  model 
size.  In  addition,  for  a  fixed  KQ,  A,  and  y,  values  of  K1  can  be  determined 
to  achieve  ACQ  -  0  for  an  interference-free  drag  simulation.  This  process  is 
shown  schematically  in  Fig.  4.  This  result  applies  to  affinely  related  model 
geometries  for  a  variety  of  fineness  ratios,  wall  heights,  and  blockage  ratios. 

2.2  CONFINED  HIGH  ASPECT  RATIO  CONFIGURATIONS 

2.2.1  Asymptotic  Representations 

The  point  of  departure  is  the  transonic  lifting  line  theory  of  Ref.  13. 
This  theoretical  model  treats  the  case  of  a  high  aspect  ratio  wing  in  a  free 
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KARMAN  GUDERLEY  SIMILARITY  PARAMETER 


K  ■  K0  +  3  K1 

H"3 


+ 


ar  *  BLOCKAGE  RATIO 


Fig.  4  Schematic  of  universal  plot  of  interference  drag 

coefficient  versus  similarity  parameter  perturbation. 

field.  Referring  to  the  wing  geometry  confined  within  the  rectangular  tunnel 
shown  in  Fig.  5,  the  free  field  wing  theory  applies  in  a  limit  in  which  if 
5  is  a  characteristic  thickness  ratio  of  a  unit  chord  wing  of  span  b,  a 
transonically  scaled  aspect  ratio  parameter,  B  =  5  b  -►  ».  This  scaling 
Is  related  to  the  transverse  extent  of  the  wave  system  which  is  0(6  in 
a  Karman-Guderley  (KG)  limit  involving  6  -*■  0  to  be  described.  In  what  follows, 
if  the  width  of  the  tunnel  is  2h,  limits  are  considered  involving  B  and 
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Karman-Guderley  Limit 

The  asymptotic  expansion  for  the  velocity  potential  $  is 
jj  =  x  +  <52/3<Kx,?,2;K,A,H,B)  +  •••  (48a) 

which  is  valid  in  the  KG  limit 

1  -  M2 

x,y  =  6^3y  ,  z  =  6^3z  »  K  =  —pjy  ,  A  -  ^  fixed  as  6  -*•  0  .  (48b) 

Substitution  of  this  expansion  into  the  exact  equations  and  boundary 
conditions  leads  to  the  KG  problem  of  transonic  lifting  surface  theory 
in  which  the  equation  of  motion 

[k  -  (y+U'I’J^x  +  +  <f>22  =  0  (49a) 

is  to  be  solved  subject  to  the  wing  and  wall  boundary  conditions 


£  U.0...Z)  ■ (x<f )-A 

(49b) 

(x,y,±H)  =  (x,±$H,z)  =  0  , 

3y  3Z 

(49c) 

where  the  equation  of  the  wing  is 

y  =  6Fu  ^(x,z/B)  -  A  ,  -B  <  z  <  B  . 

In  addition  to  (49b),  the  KG  formulation  also  involves  wake  and  far  field 
conditions  which  will  not  be  considered  in  this  discussion. 
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Inner  Region 

For  the  inner  or  near  field  region,  the  flow  is  almost  two  dimensional, 
with  corrections  due  to  the  finite  aspect  ratio  and  wind-tunnel  walls.  The 
asymptotic  expansions  expressing  this  structure  are 

4>(x»y,z)  =  <|>0(x,y;z*)  +  ^  ^  (x,y;z+)  +  (50a) 

K  =  Kq  +  f  K,  +  —  (50b) 

A  *  A0  +  E  A1  +  ***  *  (5°c) 

The  expansion  (50a)  applies  in  the  limit 
z  H 

z*  =  |,x,y,^-  =  y  fixed  as  5  -*■  0,  B  -*■  ®  independently  .  (50d) 


The  equations  for  the  approximate  quantities  associated  with  (50)  are: 

I  Kq  -  (y+1  )4>0  1*0  +  4>q_  ■  0 

L  xJ  xx  yy 

K)  -  (y+1)4>q  |4>i  -  (y+D4»i  <l>o  +  <j>i  “  - 

L  xJ  xx  x  xx  yy  xx 


which  are  to  be  solved  subject  to  the  boundary  conditions 


(51a) 

(51b) 


3<f>0 

8y 


y=0 


34'1 


y=0 


(x,z/B)  -  Aq  , 


(52a) 


(52b) 
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It  Is  to  be  noted  here  that  the  essential  interference  problem  involves 
the  solution  of  (51b)  which  is  the  wall  and  aspect  ratio  correction  to  the 
(zero  subscript)  unconfined  two-dimensional  near  field. 

To  complete  the  formulation  for  the  perturbations,  a  suitable  far 
field  is  required.  This  comes  from  matching  with  the  outer  solution.  To 
dominant  order,  this  represents  a  bound  line  vortex  shedding  its  vortlcity 
downstream  as  an  idealization  of  the  wing.  If  K  is  sufficiently  large, 
the  far  field  is  subsonic  and  is  described  by  a  three-dimensional  form 
of  the  Prandtl-Glauert  equation.  In  what  follows,  only  this  situation 
will  be  considered.  As  contrasted  to  the  free-field  case  discussed  in 
Ref.  13,  the  vortex  sheet  is  reflected  In  the  walls. 

Outer  Expansion 

The  asymptotic  expansion  exhibiting  the  anticipated  far  field  flow 
features  previously  discussed  is 

4>  •  ^g(x*,y*,z*)  +  ^ (x*,y* ,z*)  +  g-  <P2(x*,y*,z*)  +  •••  ,  (53a) 


which  Is  valid  in  the  limit 

x*  -  -  ,  y*  =  £  ,  z*  =  jr  fixed  as  B  +  “  . 

B  B  B 

The  representation  (53a)  leads  to  the  hierarchy  of  equations. 


+  +  'Pn  ^ 

U  «.*,*  U7*T* 


x*x*  y*y*  z*z 


**1**1—**.  -° 


x*x*  y*y*  'z*z* 


0*2  . . +  .  *  *2 .  .  =  *att-  Kl*0 

x*x*  y*y*  z*z*  x*  x*x* 


x*x* 


(53b) 


(54a) 

(54b) 

(54c) 
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which  are  subject  to  the  wall  conditions, 


5<4),1,2 


3y* 


^0,1,2 


y*=±3u 


32* 


Z*=±y 


=  0 


11  =  H/B 


(55) 


and  appropriate  boundedness  conditions  at  x*  =  «,  as  well  as  matching  with 
the  far  field  of  the  inner  solution  as  x*,y*  ■+■  0,  In  what  follows,  results 
from  our  first  efforts  to  obtain  this  behavior  are  indicated. 


Investigations  of  Inner  Limit  of  Outer  Flow 

Using  Green's  theorem,  the  flow  over  the  wing  and  its  wake  can  be 
shown  to  be  the  superposition  of  a  surface  source,  and  doublet  distribution 
as  well  as  a  nonlinear  volume  source.  The  doublet  suitable  for  this 
representation  satisfies  the  homogeneous  Neumann  conditions  appropriate 
for  the  solid  walls.  If  G  Is  the  potential  of  this  doublet,  dropping 
star  superscripts,  and  introducing  the  variable  X  =  x/./k^, 

AG  *  5'(y)6(z-c)6(X)  ,  (56a) 


with 


SG 

ay 


y=6u=v 


=  o  . 


Introducing  the  Fourier  transform  of  G,  i.e., 
5  =  -p=f  elkXG(X,y,z)dx  , 

(56a)  implies 

k* G  =-Lr61(y)6(z-0  . 

3y  \[Zt( 


(56b) 


(57) 
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2,2.2  Free  Field  Case 

Suppressing  the  Neumann  conditions,  (56b),  the  question  of  whether 
the  formulation  (56a)  leads  to  the  free  field  classical  Prandtl  lifting 
line  formulas  has  been  answered.  This  is  important  in  establishing  a 
“datum"  term  in  obtaining  a  correction  to  the  free  field.  Denoting  G  =  T 
for  the  free  field  solution,  the  appropriate  solution  satisfying  (57)  with 
the  proper  far  field  decay  is 


T  =  (2tt) 


-3/2  _3_ 

sy 


-  K0{kr')| 


kK£(kr')y 
/  y2  +  (z-c)2 


where  r“ 


'/y2  +  (z-?)2.  From  Inversion  and  Green's  theorem,  noting  that 


ikXKg(kr)dk  = 


9 


the  potential,  ^  q  of  a  bound  vortex  shedding  a  sheet  downstream  with  a 
shedding  strength  -y(c)  is  given  by  the  relationship 


^°»°  =  /  ^  /0  d^)T(x,y,z;S,5)  , 


where 


T(x,y,z) 


1 

(2*) 2 


kK^(kr")y 
%/  y2  +  (z-?)2 


dk  . 


(58) 


(59) 


Thus, 


X2  +  y2  +  (z-;)i 


y _ 

y2+  (z-c> 


(60) 
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which  agrees  with  the  Biot-Savart  low  speed  expression  for  this  vortex 
assemblage. 


2.2.3  Confined  Case 

For  finite  H,  the  appropriate  expression  for  G  solving  the  Sturm  Liouvllle 
problem  (57), (56b)  is 


G 


y> 

v^6„ 


coshBn(s-p)coshBn(z+p) 
's'fnhB  u 


sinXny 


where 


B_  =  A*  +  k* 
n  n 


(61) 


(62) 

(63) 
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3.0  CONCLUSIONS 

Asymptotic  theories  have  been  Investigated  for  the  treatment  of 
transonic  tunnel  wall  interference.  Two  limiting  cases  have  been  considered. 
The  first  which  was  emphasized  in  this  study  involves  low  aspect  ratio 
shapes  where  the  characteristic  lateral  dimension  of  the  model  is  small 
compared  to  its  length,  and  the  tunnel  height  is  inversely  proportional 
to  this  lateral  dimension.  Under  these  circumstances,  key  findings  are: 

1.  The  wall  correction  for  lift  is  of  higher  order  than  the  drag. 

This  contrasts  to  the  two-dimensional  case  described  in  Ref.  10 
of  an  airfoil  between  solid  walls  for  which  both  corrections  are 
of  the  same  order. 

2.  A  similarity  law  has  been  derived  for  the  correction  of  the  drag 
coefficient.  It  Implies  that  the  change  in  this  quantity  (referred 
to  the  maximum  cross  sectional  area)  due  to  solid  cylindrical  wall 
interference  is  proportional  to  the  quotient  of  the  blockage  ratio 
and  the  normalized  wall  height  H. 

3.  The  theory  which  is  applied  to  the  case  where  the  model  far  field 
is  subsonic  gives  explicit  results  for  the  rate  of  decay  of  the 
model  and  wall  perturbations  of  the  free  stream  which  are  not 
readily  accessible  from  purely  computational  correction/assessment 
simulations.  Also,  it  allows  extrapolations  to  zero  model  size 

to  be  made, 

4.  On  the  basis  of  the  theory,  a  perturbation  in  the  Mach  number  and/or 
similarity  parameter  appearing  in  the  theory  can  be  determined  to 
eliminate  the  tunnel  wall  perturbation  to  the  drag.  As  in  the  lift 
correction  for  the  two-dimensional  case,  a  universal  curve  of  drag 
interference-free  perturbations  to  the  free  stream  similarity 
parameters  as  a  function  of  the  latter  can  be  plotted  based  on  Item  2. 
This  curve,  which  is  a  generalization  of  the  transonic  area  rule, 
applies  to  aircraft  of  differing  cross  sectional  shape  variations 
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in  the  streamwise  direction  but  with  the  same  longitudinal  area 
progression.  The  blockage  ratio  and  tunnel  height  are  separated 
out  in  this  universal  variation.  The  result  of  universality  is 
that  it  achieves  major  computer  savings  in  assessment/correction 
evaluations. 

5.  The  interference  pressure  on  the  body  depends  only  on  the  streamwise 
coordinate  to  dominant  order. 

The  second  case  considered  involves  large  aspect  ratio  confined 
wings.  If  6  is  the  wing  thickness,  b  is  a  span  dimension,  and  h  is  the 
tunnel  height,  both  in  units  of  the  chord,  a  limit  within  transonic  small 
disturbance  theory  was  studied  for  which  B  =  b6^3,  H  =  h6^3  -*•  <*>.  This 
process  leads  to  a  generalization  of  lifting  line  theory  considered  in 
Ref.  13.  The  structure  of  the  flow  near  the  wing  (inner  solution)  and 
that  away  from  it  (outer  solution)  has  been  characterized.  Near  the  wing, 
it  retains  the  two-dimensional  character  of  the  free-field  case.  In  the 
far  field,  it  is  the  reflection  in  the  walls  of  a  vortex  assemblage 
consisting  of  a  bound  vortex  perpendicular  to  the  flow  shedding  trailing 
vortices.  The  intensity  of  the  latter  is  proportional  to  the  spanwise 
loading.  The  assemblage  induces  incidence  corrections  to  the  two-dimensional 
near  field.  A  part  of  this  effect  is  due  to  wall  corrections  not  present 
in  the  free-field  case. 
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It  is  recommended  that  the  theory  for  slender  bodies  be  extended  in 
several  ways.  To  treat  magnetic  support  systems,  Eq.  (9a)  can  be  suitably 
modified  by  discarding  the  dominant  term  Aq/R.  This  change  gives  the 
appropriate  far  field  doublet  rather  than“source  behavior.  Equations  { 1 2d) 
and  (12e)  by  appropriate  multiplication  by  k  give  associated  differentiation 
of  the  reflected  source  to  obtain  the  doublet.  The  same  process  previously 
discussed  will  furnish  the  asymptotic  approximation  of  these  integral 
representations  near  the  origin,  providing  the  far  field  in  Problem  PI. 

To  obtain  the  correction  for  pressure  boundary  conditions  and  other  wall 
simulations,  suitable  Green's  functions  have  to  be  employed  in  the  Green's 
theorem  leading  to  far  field  solutions  such  as  (9a).  Correspondingly, 
analogous  results  replacing  Eqs.  (12)  can  be  obtained  using  transform 
methods  and  a  convolution  theorem.  For  the  matching,  all  of  the  Ideas 
presented  herein  for  the  solid  wall,  sting  supported  case  will  be  applicable. 
In  the  case  of  rectangular  walls,  if  8  is  fixed,  it  is  anticipated  that 
the  reflections  providing  the  far  field  in  an  analogous  problem  for  PI 
will  involve  elliptic  functions. 

The  large  aspect  ratio  analysis  should  be  completed  by  determining 
the  inner  limit  of  the  outer  solution.  To  obtain  this  behavior,  an 
integral  representation  is  required  for  the  series  (61).  One  relevant 
method  considers  this  sum  as  a  residue  expansion  of  a  contour  integral. 

From  this  development,  the  procedures  described  herein  for  the  low  aspect 
ratio  case  can  be  applied.  For  the  latter,  quantitative  results  should 
be  obtained  from  a  computational  solution  of  the  problem  formulated  in 
Eqs.  (6c),  (47),  (20b),  (13c‘),  and  (9b).  Other  generalizations  besides 
the  ones  described  previously  in  this  section  involve  the  choked  case, 
and  the  K  ■+  0  limit.  These  are  of  great  practical  interest  and  should 
also  receive  attention.  The  asymptotic  methods  employed  in  this  report 
can  be  advantageously  used  to  treat  these  situations. 

Once  these  areas  have  received  attention,  the  combined  computational 
and  asymptotic  theories  should  be  applied  to  adaptive  walls  and  incorporation 
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of  viscous  effects.  Another  interesting  and  important  aspect  concerns 
consideration  of  simplifications  of  the  two  variable  methods  of  Refs.  4 
and  5  using  the  tools  developed  in  this  effort. 

On  the  basis  of  the  foregoing  discussion,  the  following  is  a  summary 
list  of  recommended  future  follow-on  areas  to  the  effort  described  in 
this  report: 

1.  Obtain  computational  solution  for  transonic  slender  configurations 
in  circular  cross  section  solid  wall  tunnels  and  apply  to  wind 
tunnel  wall  interference  assessment  and  correction  procedures. 

2.  Complete  high  aspect  ratio  analysis  and  develop  code. 

3.  Extend  analyses  to  porous,  slotted,  and  pressure  boundary  conditions. 

4.  Extend  analyses  to  other  wall  cross  sections. 

5.  Treat  choked  case. 

6.  Apply  theory  to  adaptive  walls. 

7.  Treat  viscous  effects. 

8.  Exploit  asymptotic  and  other  theoretical  methodology  to  simplify 
the  two-variable  method  of  Refs.  4  and  5. 
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APPENDIX  A 

FAR  FIELD  CALCULATION  OF  $0 

To  determine  the  behavior  of  4>q  for  large  R,  Green's  theorem  Is 
applied  to  the  region  shown  in  Fig.  A1 .  For  this  purpose,  (6a)  is  written 
in  a  "stretched"  coordinate  system  involving  the  variable  X  defined  in 
(8b).  This  gives  Poisson's  equation, 

«o  ■  Ko3/2<^'>yoxx  ■ 


where 


Fig.  A1  Region  for  application  of  Green's  theorem. 
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In  addition,  a  fundamental  solution  G  is  defined  with  the  property  that 
AG  =  SCx-S)  6(0-8-)  .  (AZ) 

Green's  theorem  applied  to  the  region  shown  in  Fig.  A1  can  be  written  as 


(A3) 


where  n  denotes  differentiation  in  the  direction  of  the  outward  normal  to 
the  region  V-j .  In  (A3),  the  left  hand  side  represents  an  integration  over 
the  surfaces  S  ,  S  ,  and  S^.  The  surface  S  is  a  slender  cylinder  of 
radius  r  =  e  with  the  X  axis  as  its  axis  of  symmetry.  The  surface  is 
the  total  area  of  the  shocks,  and  Sm  is  an  infinitely  large  sphere. 

The  right  hand  side  of  (A3)  represents  a  volume  integration  over  V-| . 
On  the  basis  of  (A2), 


G  = 


_ l _ 

'J  (x-s)2  +  r2  -  2rp  cos(e-e') +  p2 


(A4a) 


3G 

ap 


1 

4? 


p  -  r  cos(9-8 1 ) 


[(X-€)2  +  r2  -  2rp  cos(9-6 ■ )  +  p2] 


37? 


(A4b) 


For  a  sting  support,  it  is  anticipated  that  will  behave  dominantly  in 
a  subsonic  far  field  like  an  unconfined  source,  l.e.,  <J>Q  -  -  .  From 

(A3)  and  this  assumption,  the  contribution  from  is  0{R-1)  as  R  +  “ 
and  therefore  vanishes.  For  the  magnetic  suspension  case,  the  doublet 
behavior,  <J>q  =  0(R  2)  is  expected  and  again  the  contribution  from  SM 
vanishes.  The  integral  on  the  shock  surfaces  also  disappears.  This 
can  be  demonstrated  in  the  same  manner  as  in  the  two-dimensional  case  by 
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Integration  by  parts  of  the  right  hand  side  of  (A3).  Combination  of 
boundary  terms  from  the  partial  integration  and  the  So  integral  nullifies 
the  shock  contribution  by  the  shock  relations. 

Implementing  this  procedure,  (A3)  can  be  written  as 


h  -  r3  +  '4  +  'a 


where 


I 


1 


I 


2 


I 


3 


.  (A5) 


(A6a) 


(A6b) 

(A6c) 


(A6d) 


(A6e) 


On  the  basis  of  anticipated  matching  with  the  Axis  Layer  as  described  in 
connection  with  (24c), 


S'  (x) 


2rr 


Jin  r  +  g(x) 


as  r  -*•  0  . 


(A7) 
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JVp 

Approximating  G  and  for 


0,  and  use  of  (A7)  gives 


I 


1 


&ne  +  g(£) 


e(-r  cos (8-9 1 )  +  ■»«)d9‘ 
4ir[(X-?)2  +  r2-0(e)]3/2 


=  0(e  £ne)  as  e  ■+  0  . 

In  a  similar  manner,  noting  that  S'(x)  =  0  for  x  >  1, 


(AS) 


*2  “  ‘  4?  J  TT"1"  =T=  ‘  W  ^  J1  * 

*  ■'o  y  R2  -  2£R  cosw  +  £2  q  n=0  Rn  1  •'O 

from  expanding  the  denominator  in  Legendre  polynomials  and  termwise  integration. 
On  integration  by  parts  and  letting  R  +  »,  this  becomes 


I  =  SH)  + 
lZ  W 


fs(  1 }  -  f1  S{e)d€j£2^+(s(l)-2  f1  5S(5)dJ 

[  •/0  )  4ttR2  (  J0  ) 


PgCcOSto) 

4uR3 


+  CJ(R-1*) 


(A9) 

where  the  factor  in  the  second  braces  representing  the  strength  of  the 
quadrupole  contributes  to  the  quantity  A  in  (9a). 

By  (A2),  (A6c)  gives 

I3  =  V-  (A10) 

From  (Al),  (A6d)  becomes 

U  -  •  (All) 

Integrating  (All)  by  parts,  combining  the  boundary  terms  with  those  in  1^, 
applying  the  shock  relations,  and  collecting  the  results  for  1^,  Ig,  and 
Ig  gives  the  first  two  terms  of  (9a).  These  are  the  first  two  spherical 
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harmonic  solutions  of  (Al),  assuming  the  response  to  the  nonlinear 
forcing  term  is  negligible.  The  higher  order  terms  include  the  response 
to  this  forcing  term.  For  this  purpose,  the  forcing  term  is  evaluated 
for  large  R,  from  the  dominant  representation  of  <fg. 

Accordingly, 


2A, 


0 

<fin  = - -  P1(cosw)P,(cosa>) 

•  nS  I  £ 


A080 


'X  XX 


(6P1P3  +  P|)  + 


( Al  2 ) 


where  from  (9a), 

Aq  =  A0//T0  ,  Bg  =  Bq/Kq  , 

and  Pn(costu)  is  a  Legendre  polynomial  of  order  n.  The  products  of  the 
Legendre  polynomials  can  be  expressed  as  a  sum  of  them  (Ref.  Al).  This 
fact  facilitates  the  determination  of  the  particular  solutions.  In 
particular,  denoting  u  =  cosoj. 


P-,(u)P2(u) 


3P3  +  2P1 
5 


(A13) 


To  obtain  the  particular  solutions,  (Al)  is  written  in  spherical  coordinates, 
and  the  method  of  separation  of  variables  is  used.  Considering  the  response 
to  a  typical  forcing  term  in  (A12),  using  (A13),  this  gives 


-  1  3 


A$0  a 


&R 


Pn(cOSd)) 

R5 


fail 

Gradshteyn,  I.S.  and  Ryzhik,  I.M.,  Tables  of  Integrals,  Series  and 
Products,  Corrected  and  Enlarged  Edition,  Academic  Press,  New  York, 
Eq.  (Byib,i>)’  p.  102b. 
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In  which 


♦0  =  *(R)T(w)  . 


(A14) 


If  T  =  Pn(cosco} , 


R~a 

6^“n(n+l)  ’ 


(ATS) 


The  special  case,  n=2  in  (A15)  leads  to  irrelevant  logarithmic  solutions. 
Others  are  of  higher  order  and  are  indicated  by  the  0  symbol  in  (9a).  These 
are  unimportant  for  the  matching  discussion  In  this  report.  On  the  basis 
of  this  analysis,  the  forcing  term  response  particular  solution  to  the 
dominant  source  term  of  <t>0»4>g  is 


- |cos3(d  +  l  cosul  , 

16ff2K0[x2  +  K0?2]3/2  '  5  ' 


(A16) 


and  leads  to  (9a). 

In  connection  with  this  derivation,  it  should  be  noted  that  in  (A8) 
the  lift  effect  is  associated  with  a  doublet  term  (not  shown),  n*  cose'/e. 
that  is  in  the  bracketed  factor  in  the  integrand.  This  term  is  higher 
order  In  5  in  the  limit  (4)  as  can  be  inferred  from  matching  of  the  Axis 
and  Central  Layers.  Therefore,  it  has  been  neglected  in  this  analysis. 
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APPENDIX  B 

DERIVATION  OF  DOMINANT  APPROXIMATION  FOR 
WALL  LAYER  <PQ  GIVEN  BY  EQ.  (12a) 


The  singular  behavior 


as  R+  0  , 


R+ 


R/H 


9 


(Bl) 


by  virtue  of  the  need  to  match  with  (j>Q,  is  anticipated.  Accordingly*  (11a) 
can  be  rewritten  as 


u  /kT  2itrT 


(B2) 


where  I  6  (r)dr  =  1.  The  appropriate  boundary  condition  associated  with 
•'O 

solid  walls  is 


(B3) 


Applying  the  divergence  theorem  to  the  region  bounded  by  a  cylinder 
consisting  of  the  curved  surfaces,  r^=l,  and  the  flat  faces  X  ■*  «,  (B2) 
and  (B3)  imply  that 


_  S(l)  X 


(84) 
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To  solve  the  boundary  value  problem  embodied  in  (B2)  and  (83),  -the 
exponential  Fourier  transform  is  used.  For  the  purpose  of  obtaining  . 
convergent  integrals,  the  property  (B4)  is  utilized  to  define  a  regularized 
version  of  v>q  denoted  by  M  for  which 


X+ 


(B5) 


Accordingly, 


AM 


UV)  1 

I  2irr*  * 


a(x)  . 


(B6) 


With  the  following  exponential  transform  pair. 


(B7) 


(B8)  ' 


the  subsidiary  equations  for  M  become: 


1  d_  L  dM \  .  .  6+(r)  _  I 

r  dr  \r  Hr/  2nr  t  * 


(B9) 


The  first  forcing  term  on  dropping  the  daggers  on  r  and  X  can  be 
eliminated  and  replaced  by  the  boundary  condition 


11m 
r-*-  0 


dM  = 
dr  ~  2 rr 


(BIO) 
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The  other  boundary  condition  is 


dM 

3r 


=  0 


r=l 


(Bll) 


The  solution  of  the  boundary  value  problem  is 


'K^(k) IQ( kr)  -  Kq{ kr) I^( k) 


"2TTt‘(k) 


dk 


uka 


(B12) 


where  the  effect  of  the  last  term  regularizer  in  the  braces  eliminates  the 
double  pole  at  the  origin.  Since  Iq(z)  =  Ij(z)  has  only  pure  imaginary  simple 
zeroes,  the  appropriate  inversion  contour  is  a  large  semicircle  with  its  base 
along  the  real  axis.  For  |k[  +  the  square  bracket  term  Is 
e-|u|r-TX 


,  where  k  =  p+iT.  Accordingly,  the  semicircle  is  in  the  upper 


half  plane  for  X  >  0,  and  in  the  lower  for  X  <  0.  The  zeroes  of  Iq  are  at 
k*±iXn,  n=l,Z,*«».  is  the  solution  of  the  equation 

J-, (An)  =  0  .  (B13) 

On  summing  the  residues  at  the  poles  corresponding  to  these  zeroes,  we  obtain 


XI 


JQUnr) 


Cd0(An)]2 


( B1 4) 


where  the  Wronskian  relation 
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and  (B13)  have  been  utilized  in  (B14),  which  can  also  be  derived  by 
eigenfunction  expansions. 

Equations  (12c)-(12e)  follow  from  (B12)  and  the  properties  of  the 
modified  Bessel  functions*  upon  restoring  the  daggers  and  noting  that 
(B9)  admits  solutions  for  positive  and  negative  values  of  k. 
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APPENDIX  C 

DERIVATION  OF  BEHAVIOR  OF  <PQ  NEAR  ORIGIN  R+=0  GIVEN  BY  EQ.  (13a) 

To  obtain  the  required  representation,  { 1 2d)  and  (12e)  will  be 
considered.  Now,  Mq  can  be  evaluated  exactly  from  a  table  of  Fourier 
transforms  to  give 


M 


0 


1 

4*Rf  ' 


(Cl) 


This  Is  precisely  the  source  behavior  desired  for  matching. 

One  candidate  method  of  evaluating  Mj  involves  expansion  of  the 
Integrand  near  the  origin.  This  leads  to  divergent  Integrals.  To  avoid 
this  difficulty,  the  derivatives  of  M1  with  respect  to  rf  and  X+  are 
considered  using  the  same  procedure.  Accordingly,  on  dropping  the  daggers. 


-/0”i2^Ki(k)Ii(kr>‘,k 


«  -  rl 


1 


> 


(C2a) 


on  expanding  the  trigonometric  factor  cos  kX  and  integrating  term  by 
term  with 


k2Uk) 

T^rdk  ■ 


(C2b) 


_  -2k 

Since  the  integrand  of  I-j  is  bounded  at  the  origin  and  0(k2e"  )  as  k  +  ®, 

it  is  evident  that  1^  converges. 

Moreover, 
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— K,  { k) 


i 


(C3) 


Again,  expanding  the  trigonometric  factor  and  recognizing  that  the  last 

integral  Is  well  known,  and  can  be  evaluated  by  Cauchy's  theorem  as 
tt  3M, 

—  sgn  X,  -ipg-  is  given  as 


3M1  X  -  TT 

^  f  I,  -  I  sgn  x  .  (C4) 


Since 


which  is  itself  a  convergent  integral,  utilization  of  (C2),  (C4),  and  (C5) 
with  (12a)  leads  to  (13).  This  result  could  also  have  been  obtained  from 
a  Taylor's  expansion  of  M-j  near  the  origin.  It  is  interesting  to  note 
in  this  connection  that  the  second  term  In  (C4)  leads  to  cancellation  with 
the  first  term  in  (12a)  as  -*■  0. 
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A 

A 

Ag,i 

A_0 

a0 

ao 

ar 

B 

B_0 

b 

b0 

bo 

lo 

CD 

AC0 

CP 


D 


F 


G 


g 

q*(n) 

90 

g„(x) 


•  NOMENCLATURE 

Angle  of  attack  parameter 
Constant  appearing  in  (9a) 

Approximations  in  angle  of  attack  parameter  expansion  (5d) 
Constant  appearing  in  (9a) 

Constant  appearing  in  (13b) 

Constant  defined  in  (22b) 

Blockage  ratio 

1  /3 

Reduced  span  =  5  '  b,  body  function 
Constant  appearing  in  (9a) 

Span 

Constant  defined  in  (13c) 

Constant  defined  after  Eq.  (20b) 

Constant  appearing  in  (9a) 

Drag  coefficient 

Wall  interference  effect  on  drag  coefficient 

Pressure  coefficient 

Chord 

Drag  force 

Body  shape  function  in  (1) 

Green's  function 
Function  defined  in  (44) 

Functions  appearing  in  (26a) 

Functions  appearing  in  (30) 
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H 

h 

I. 


Kq.K] 


V1! 

Pn(costu} 

q 

R 

R 

R+ 

R 


r 


r* 

S(x) 

S_ 


Reduced  tunnel  radius  in  (4),  and  tunnel  height  defined  before  (48a) 
Tunnel  dimension  in  units  of  chord 

Modified  Bessel  function  and  designation  of  integrals  in  (A5) 
Integral  In  (A5) 

Transonic  similarity  parameter  in  Eq.  (4) 

Modified  Bessel  function 

Dominant  two  terms  in  asymptotic  expansion  for  K  (5c) 

Fourier  transform  variable  of  integration 
Portions  of  ^  solution  in  (12a) 

Legendre  polynomial 

Dynamic  pressure 

Reduced  polar  radius  in  (8b) 

Scaled  polar  radius  in  (8a) 

Reduced  polar  radius  =  R/H 
Intermediate  variable  R/n(H) 

Intermediate  variable  R/c(H) 

Function  of  R  in  separation  of  variables  solution  (A15) 

Cylindrical  coordinate 

5r 

r/H 

r/5 

Reduced  cross  sectional  area 
Area  in  Fig.  A1 
Area  on  shock  in  Fig-  A1 
Area  at  infinity  in  Fig.  A1 
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&>{*)  Source  strength  =  S'(x)/2tt 

T(u)  Function  of  co  in  separation  of  variables  solution  (AT 5) 

U  Freestream  velocity 

u  cosw 

Volume  of  region  in  Fig,  A1 

X  x//K^ 

X+  X/H 

x,y,z  Dimensional  Cartesian  coordinates 

x, y,z  Nondimensional  Cartesian  coordinates 

y, z  6y,  fiz,  respectively 

■j"  mj 

y  ,Z  y/H,  zt H,  respectively 

x*  x/B 

y*,z*  y/6,  z/S  or  y/B,  z/B 

a  Angle  of  attack 

3  Aspect  ratio  of  rectangular  tunnel  cross  section  defined  in  Fig.  5 

\  +  k.2 

e  Radius  of  internal  cylindrical  boundary  S  in  Fig.  A1 

Y  Specific  heat  ratio 

<5  Thickness  ratio,  characteristic  flow  deflection 

<5(x)  Delta  function 

A  Laplaclan  operator 

en  Gauge  function  appearing  in  Wall  Layer  expansion  {10) 

S,n,C  Duirnny  Cartesian  coordinates 

0  Azimuth  angle  in  spherical  coordinates,  cylindrical  coordinates 

0'  Dummy  variable  for  0 
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«n  Gauge  function  appearing  in  angle  of  attack  parameter  expansion  (5d) 

Ap  Eigenvalue  appearing  in  (63)  and  (B13) 

y  H/B  in  (55) 

Gauge  functions  appearing  in  Central  Layer  expansion  (5a) 
vn  Gauge  functions  appearing  in  similarity  parameter  expansion  (5c) 

v  Bp 

w  Polar  angle  in'  spherical  polar  coordinates 

p  Dummy  variable  for  r 

yfi  Approximations  for  perturbation  potential  appearing  in  Wall  Layer 

n  expansion  (10)  and  (53a) 

<J>  Perturbation  potential  defined  in  (48a) 

cfi  Perturbation  potential  approximations  in  Central  Layer  representation 

n  (5a)  and  (50a) 

<j>*  Perturbation  potential  approximations  in  Axis  Layer  representation  (23a) 

$  Velocity  potential 

Tn  Gauge  function  appearing  in  (23a) 

Subscripts 
B  On  body 

H  Lower  surface 

u  Upper  surface 

Special  Symbols 

Tilde*1"  }  Denote  Fourier  transform 
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